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We present results from ah initio calculations of the mechanical properties of the rhombohedral 
phase (/3) of vanadium metal reported in recent experiments, and other predicted high-pressure 
phases (7 and bcc), focusing on properties relevant to dynamic experiments. We find that the volume 
change associated with these transitions is small: no more than 0.15% (for /3 - 7). Calculations of 
the single crystal and polycrystal elastic moduli (stress-strain coefficients) reveal a remarkably small 
discontinuity in the shear modulus and other elastic properties across the phase transitions even at 
zero temperature where the transitions are first order. 

PACS numbers: 62.50.-(-p, 61.50.Ks, 31.15.Ar, 62.20.Dc 
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I. INTRODUCTION 



The existence of a high-pressure rhombohedral phase 
of pure crystalline vanadium has been the focus of an 
intense research effort recently. The first indication of 
a phase transition came from the theoretical observa- 
tion that the C44 shear modulus of bcc vanadium de- 
creases and becomes negative at pressures greater than 
~1.3 Mbarjii^ pressures that are experimentally acces- 
sible. A negative shear modulus means the material is 
mechanically unstable under trigonal (prismatic) shear, 
suggesting a phase transition. At that time the ex- 
perimental evidence showed no phase transition up to 
1.54 Mbar.'^ Then recently, Mao and coworkers^ con- 
ducted x-ray diffraction experiments in the diamond anvil 
cell (DAC) up to 1.5 Mbar and found features in the 
diffraction peaks that were consistent with a second- 
order phase transformation to a rhombohedral struc- 
ture with an R3m point group symmetry at pressures 
above 0.69 Mbar. It was soon confirmed that density 
functional theory (DFT) finds the rhombohedral phase 
to be the ground state at zero temperature and pres- 
sures above 0.8 Mbar, in reasonably good agreement with 
experimenti^ In fact, it was shown that DFT predicts ad- 
ditional phase transformations that had not been found 
in experiment, i.e. a first-order transformation to a dif- 
ferent rhombohedral structure at 1.2 Mbar and a third 
transformation back to the bcc structure at 2.8 Mbari^ 
The prediction of the existence of the two high pressure 
phase transformations has been subsequently confirmed 
with DFT phonon calculations 1^ 

Alternative techniques may provide the pressures 
needed to observe the second rhombohedral phase and 
the reentrant bcc phase. Dynamic experiments do not 
rely on the mechanical integrity of anvils and are able 
to reach multi-megabar pressures. They have been used 
to study similar transformations such as the diffusion- 
less a-e transition in iron.^ There are several challenges 
specific to vanadium, however. The softening of the 
shear modulus and the rhombohedral phase transition 
are related to subtle electronic effectsr' which are likely 



to be weakened by increased temperature. Recent ramp 
wave techniques based on laseri^iiS and graded-density 
impactor— drives are able to generate high pressure with- 
out the entropy generation of Shockwave techniques, and 
are therefore preferable in the present context. An- 
other challenge is that the subtle rhombohedral distor- 
tion (< 1°) detected by x-ray diffraction in the DAC is 
probably too small for in-situ x-ray diffraction in dynamic 
experiments.^^ Indirect techniques are an alternative to 
detect the transition. For example, VISAR free-surface 
velocity measurements can detect changes in the density 
due to a volume change, and they can be used to infer 
the longitudinal stress, and hence the change in strength 
if the equation of state is known i^°d^ Rayleigh- Taylor 
growth rate is another way to probe strength.^ 

In this article we use DFT to make predictions about 
the properties of high-pressure vanadium relevant to dy- 
namic experiments. We compute the magnitude of the 
volume change associated with the three phase tran- 
sitions related to the rhombohedral structure in Sec- 
tion |TTT1 We also compute the elastic properties and cal- 
culate bounds on, and an explicit estimate of, the poly- 
crystalline shear modulus in Section |TV] and Section |V] 
respectively. Since the strength is typically assumed to 
vary with the shear modulus^ any anomalies in the 
shear modulus are likely to provide a signature in the 
VISAR trace. Indeed, an important motivation for the 
present work is to assess whether the bcc shear modulus 
C44 going to zero is likely to produce a strong signa- 
ture. The shear modulus also affects defect energetics, 
and may have a measurable effect on transition kinetics. 
We consider the implications of our results for dynamic 
experiments to detect the high-pressure phases. 



II. THEORETICAL BACKGROUND 

The rhombohedral crystal structure of vanadium at 
high pressure results from a slight distortion of the bcc 
structure. Specifically the distortion is a uniaxial strain 
along (111), which remains a three- fold symmetry axis of 
the crystal. This structure is known as the /3-Po structure 
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[Strukturhericht Ai, Pearson hRl). It still has a single 
atom per unit cell, so the rhombohedral transition may 
be expected to be difFusionless (martensitic) and likely 
rapid despite the small energy difference. There are four 
independent three- fold axes, so there are four variants of 
the rhombohedral crystal that are degenerate in energy. 

The ground state of the single-crystal rhombohedral 
phase has been determined from first principles using a 
volume-conserving rhombohedral shear pathj^ 

(k 6 5\ 
T{5) k 5] (1) 

\5 6 k j 

in the usual bcc crystal frame where k is determined from 
the real positive solution of det(T) = 1 to insure constant 
volume. The approach is to use DFT in combination 
with a gradient-corrected exchange and correlation en- 
ergy functionali^ as implemented in the Vienna Ab-initio 
Simulation Package (VASP) code along with the projec- 
tor augmented-wave (PAW) method. Specifically, the 
PAW potentials with 13 valence electrons (3s, 3p, 3(i, 
and 4s states) are used. The planewave cutoff energy is 
66.15 Ry and an unshifted 60x60x60 uniform mesh is 
used for the fc-point sampling: this results in 5216 and 
18941 k points in the irreducible Brillouin zone of the un- 
strained bcc and rhombohedral lattices, and up to 18941 
and 54932 k points for the strained bcc and rhombohe- 
dral lattices respectively. For all of the calculations we 
use a primitive cell. 

III. VOLUME CHANGE DUE TO 
TRANSFORMATION 

In Ref. H, we calculated the enthalpy and pressure 
as functions of strain along the rhombohedral deforma- 
tion path, and used the enthalpy to find any stable or 
metastable crystal structures. We noted that the equa- 
tions of state (EOS) for the bcc and rhombohedral struc- 
tures are nearly identical, so their bulk moduli are essen- 
tially equal (differing by no more than 3%), and reported 
the EOS of the ground state up to 2 Mbar. We now use 
those data together with additional data on the EOS of 
the metastable structures to calculate the volume change 
associated with the phase transformations in a readily ac- 
cessible form. 

Using the EOS Pi{V) for the stable and metastable 
structures (i = bcc, /3, 7) , we have calculated the associ- 
ated volume change AF according to 

P,{V, + ^V) ^ P,{Vi) (2) 

for pairs of structures i and j. In practice, we have calcu- 
lated the pressure at a set of volumes and used piece- wise 
quadratic interpolatation to solve the equal pressure con- 
dition ^ between those points, equivalent to the com- 
mon tangent construction at the phase boundaries. The 
EOS data are tabulated in Table [D The relative volume 



TABLE I: Equations of state for the bcc, 13 and 7 phases and 
metastable structures. Pressures are in Mbar, and volumes 
are in units of the ambient volume Vo=13.518 A''. 
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FIG. 1: (Color online) The volumetric strain AV/V of the 
ground state of single-crystal vanadium at zero temperature 
with respect to that of the bcc structure. Dashed lines cor- 
respond to three volume changes associated with three phase 
transitions. The light (dark) gray area represents the pressure 
range in which the /3 (7) phase is the ground state respectively. 
The bcc phase is stable in the other regions. 



change with respect to the bcc phase, AV/Vbcc, is plotted 
in Fig. m 

In principle, there is a volume change during the trans- 
formation from bcc (a) to the first rhombohedral phase 
(/3), and a second volume change associated with the 
transformation to the second rhombohedral phase (7), 
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and a third volume change associated with the transfor- 
mation back to bcc (a) at high pressure. There is no 
path connecting (3 and 7 that preserves the point group, 
so they must be distinct phases; the two bcc regions ap- 
pear to be connected at finite temperature.^ In practice, 
the kinetics of the transformation may cause the trans- 
formation to be overdriven so that the initial phase is 
retained in a metastable state past the phase boundary 
until the new phase has a chance to nucleate and grow. 
For this reason it is interesting to examine the entire 
curve in Fig. [U and not just the values in between the 
equilibrium phase boundaries. 

In each case the initial /SV is a volume change, so the 
volume is reduced following the transition. The volume 
change associated with the bcc to (3 transformation is 
small, 0.03% or less in magnitude. It would not be easy 
to detect such a small change in a dynamic experiment. 
The magnitude of the volume change associated with the 
second transformation is larger: about 0.15% for the (3 to 
7 transformation at 1.2 Mbar. The volume change would 
be about the same if bcc were retained to a pressure of 
~1.2 Mbar and then transformed directly to 7. This 
However, if the bcc or (3 phase persists to higher pres- 
sures, the volume change becomes progressively smaller 
and eventually changes sign, becoming a volume expan- 
sion near 2 Mbar. The final transition back to bcc again 
has a change of over 0.1% in magnitude. So the (3 to 
7 transformation has the strongest signature in terms of 
volume change, but it may not be large enough to detect. 



IV. SINGLE CRYSTAL ELASTIC MODULI 

We next consider how the single crystal elastic moduli 
change with pressure. Specifically, we calculate Bijki (-P), 
the elastic moduli with respect to a shear-stress-free ref- 
erence state at pressure P (either the bcc or rhombohe- 
dral structure, as specified). The Bij^i are often called 
the stress-strain coefficients, as we do below They are 
directly related to sound velocities at high pressure. Bijki 
can be obtained from the deformation paths used for or- 
thorhombic latticei^ or for trigonal lattice, and the de- 
tails are given in the appendix. 

The six independent stress-strain coefficients. En, 
B33, B12, Bi3, B44, and B24 (here given in Voigt nota- 
tion in the rhombohedral frame with directions [110]f,cc, 
[Il2]bcc and [lll]f,cc as 1, 2 and 3, respectively), are plot- 
ted in Fig. [IJa). They are discontinuous at the first- 
order phase transitions; however, within the domain of 
each stable phase, most of the stress-strain coefficients 
increase monotonically with pressure. The exception is 
B33 near the 7 to reentrant bcc phase boundary (roughly 
2.8 Mbar). Since ^33 is associated with uniaxial strain 
along the three-fold axis, its anomalous behavior is sug- 
gestive, but a better presentation is needed to separate 
the effects of shear and compression. We turn to it now. 

Indeed, there is a remarkable approximate continuity 
of the elastic properties across the phase transitions that 





FIG. 2: (Color online) The elastic moduli (specifically, the 
stress-strain coefficients) of the ground state single-crystal 
structure as a function of pressure; (a) Bijki written in Voigt 
notation in the frame of the primitive rhombohedral cell, and 
(b) the corresponding eigenvalues of the 9x9 stress-strain co- 
efficient matrix Bijki (see text). In the rhombohedral phase 
there are six independent elastic moduli (stress-strain coef- 
ficients) (vs. three for bcc)fi^ and three independent shear 
eigenvalues for vanadium metal (dashed curves represent dou- 
bly degenerate eigenvalues). The anisotropy ratio A"^ is also 
plotted, showing an extreme value just below the 7-bcc tran- 
sition. 



is not readily apparent from the elements of Bijki ■ The 
bulk modulus of the rhombohedral phases is within 3% 
of that of the bcc structure, as we already discussed. 
The eigenvalues of the 9x9 matrix -B(ij)(fei) provide a 
description of the elasticity that is less coordinate de- 
pendent, but there is a technical issue. In the rhombo- 
hedral phase, shear and compression are mixed in the 
sense that a non-equiaxed strain is required to produce 
purely hydrostatic pressure, and the tetragonal strains 
to produce hydrostatic pressure and pure shear are not 
orthogonal. To eliminate any ambiguity, we restrict to 
the space of constant volume strains using a projection 
matrix Il(^j){ki) = S{ij)(ki) - 5'^y4i- Then the 9 X 9 ma- 
trix ni?n/2 has 5 nontrivial eigenvalues, corresponding 
to different shear moduli. This matrix is closely related 
to von Mises stresses. 

The eigenvalues are plotted in Fig.[2Jb). The top curve 
represents two degenerate eigenvalues that are equal to 
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B' = {Bii — Bi2)/2 in the bcc phase, the usual shear 
modulus for tetragonal shear in the cubic crystal. It is 
quite smooth. The remaining three eigenvalues are de- 
generate in the bcc phase and equal to ^44, the shear 
modulus for trigonal shear in the cubic crystal (not to be 
confused with the Ba in the rhombohedral frame). In 
the rhombohedral phases two of these eigenvalues remain 
degenerate but one splits off. That single eigenvalue rep- 
resents a pure shear corresponding to the rhombohedral 
deformation. Its value is {Bn + 2^33 -I- B12 — 4:Bi^)/6 
which decreases toward zero as the pressure in the rhom- 
bohedral phases approaches the bcc phase boundary. 
This decrease is most pronounced approaching the high- 
pressure reentrant bcc phase (2.8 Mbar), but it is present 
at both. In the energy curves, it is clear that the width 
of the rhombohedral well is broadening with the change 
in pressure as it rises above the bcc well and quickly be- 
comes unstable. By the same token, the single eigenvalue 
reaches its maximum at 1.87 Mbar, the pressure that the 
rhombohedral well is deepest and most stable against the 
bcc phase. The eigenvalues can also be used to study 
the elastic anisotropy of the crystal. The anisotropy ra- 
tio {A = B44/B' in bcc) has been calculated for all four 
phases from the ratio of the smallest and largest eigenval- 
ues. For an isotropic material A = 1; for vanadium 1/A 
is never less than its ambient value, fluctuates through- 
out the entire pressure range studied, and becomes ex- 
tremely high near the 7-bcc boundary {A ^ 1/130). For 
comparison, the most anisotropic cubic transition metal 
at ambient conditions is copper^*' with A — 3.21, and 
among all cubic elements recent calculations found for 
polonium A=l/6 to 1/18 at T=OK^ 



V. POLYCRYSTALLINE SHEAR MODULI 

Polycrystalline vanadium without texture has isotropic 
mechanical behavior, described by just two independent 
elastic moduli: the bulk modulus K and the shear mod- 
ulus G. Regardless of phase, K is within 3% of that of 
the bcc structure as mentioned earlier. Using the single- 
crystal Bijki, G may be bounded by the Voigt and Reuss 
approximations of constant strain and constant stress, 
respectively. We have calculated these approximations 
using expressions equivalent to those in the literature. 
Since dynamic experiments conducted at Z-pinch and 
laser facilities often use thin-film targets with microstruc- 
tures that can vary from columnar to equiaxed depending 
on how they are processed, the Voigt and Reuss bounds 
are helpful in assessing the range of possible responses. In 
calculating the Voigt and Reuss bounds shown in Fig. [3] 
(as well as the explicit polycrystalline calculations be- 
low), we assume that the deformations are infinitessi- 
mal. With the low energy barriers, switching between 
variants of the rhombohedral phase may contribute to 
the strain with no cost in stored elastic energy, leading 
to a reduction in the shear modulus. At larger strains 
the response to rhombohedral strains stiffens anharmon- 
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FIG. 3: (Color online) The polycrystalline shear modulus 
as a function of pressure, along with the Voigt and Reuss 
bounds for bcc, /3 and 7. This shear modulus is based on the 
stress-strain coefficients Bijui of the single-crystal structure 
(see Fig. [2] (a)), and calculated from the virtual test sample 
with random grain orientation shown in color (inset). 



ically. Both of these effects have been neglected. The 
homogenized shear modulus in the rhombohedral phase 
is positive, indicating mechanical stability. The variation 
in the Voigt-Reuss difference results from the changing 
crystalline anisotropy. 

In the case of a microstructure with more equiaxed 
grains, we calculate the polycrystalline shear modulus by 
homogenizing the single-crystal B^^i of the ground state 
structure at each pressure using a virtual test sample 
(VTS).^'^ The procedure in Ref. 23 has been repeated: 
the VTS is strained in six pure shear modes (£12, £23, 
£31, £11 - £22, £22 - £33, and £33 - £11), and the results 
from six calculations are averaged to reduce the effects of 
anisotropy. The isotropic average shear modulus(Gv'Ts) 
plotted in Fig. [3| has been obtained by equating the cal- 
culated elastic energy per volume to an ideal elastic solid 
at the same strain: 

u = GvTs{e\^ + £22 + £33 + £?2 + £23 + £31) • (3) 

Only one or two of the six strains are nonzero at each run 
depending on which of the six modes has been applied, 
and the resulting variance in six independent runs is in- 
dicative of the anisotropy. The overall VTS prediction 
lies between the Voigt and Reuss bounds, and the VTS 
values are within 5% of the Voigt-Reuss-Hill average^'* ex- 
cept at the reentrant bcc phase boundary, where the VTS 
value is 31% greater. The Voigt-Reuss difference at this 
phase boundary is quite large: 1.34 Mbar. The constant- 
stress Reuss average is sensitive to the most compliant 
orientation, whereas the Voigt average is fairly insensi- 
tive; the VTS shear modulus is closer to the Voigt value. 
It significantly decreases at this point of high anisotropy, 
and may lead to an anomalous dynamic response. 
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VI. CONCLUSION 

We have investigated the properties relevant to dy- 
namic experiments for two high-pressure rhombohedral 
phases in vanadium metal. It will be challenging for dy- 
namic experiments to detect the rhombohedral phase un- 
ambiguously. The distortion is probably too small for in- 
situ x-ray diffraction, although it might be large enough 
in 71^ We have predicted that the volume change associ- 
ated with any phase transformation up to 3.15 Mbar is 
small, and may not have a clear signature in the VISAR 
trace. We have also predicted values for the single crystal 
and polycrystalline stress-strain coefficients in the rhom- 
bohedral phases at zero temperature. The (3 and 7 phases 
smoothly cut off the negative values of the bcc -B44. The 
first order transitions between the bcc and rhombohe- 
dral phases give remarkably small changes in the stress- 
strain coefficients, as evident from the plots of the shear 
modulus and the stress-strain matrix eigenvalues, apart 
from near the 7-bcc transition where the crystal is highly 
anisotropic. 

The results here were obtained using DFT at zero tem- 
perature for pure vanadium. Since the phase transition 
is driven by rather subtle electronic structure effects, the 
elastic constants may be substantially affected by changes 
in temperature or impurities j^i^i^ It would be interesting 
to see whether the remarkable continuity of the moduli 
persists as the phase boundaries and the relative stiffness 
of the bcc and rhombohedral structures change. 
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APPENDIX A: CALCULATION OF 
STRESS-STRAIN COEFFICIENTS IN THE 
RHOMBOHEDRAL LATTICE 

The high-pressure stress-strain coefficients Biju can 
be obtained in many different ways, but in Table [Til 
we summarize the deformation gradients and the corre- 
sponding strain energy relations that we used to calcu- 
late -Byfc; here. The pressure term is involved in some 
of the strain energy relations, for which the deformation 
gradients are not volume-conserving. The stress-strain 
coefficients Bijki{P) are equal to Cijki when the pressure 
vanishes, as explained in detail in Chapter 2 of Ref. [l3. 
For a recent discussion of the stress-strain coefficients 
Bijki, see Ref. [2^. For P ^ 0, the relation isM- 

Bijki = -P{SjiS^k + SiiSjk - 5ij6ki) + Cijki ■ (Al) 
The resulting stress-strain coefficients along with the cal- 
culated polycrystalline shear moduli are tabulated in Ta- 
ble |TTT1 
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TABLE II: Deformation gradients for the six independent stress-strain coefficients in the rhombohedral lattice and the corre- 
sponding strain energy relations per unit volume at pressure P. 



Stress-strain coefficient 



Deformation gradient 



Strain energy 



-B33 



B12 (and Bee) 



B24 (— S14) 



1 + S 



T{S) 



T{5) 



T{6) 
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1 











1 











1 + S 


5 








1 - (5 
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1 + 5 



\ 



i_ 

i-i 

\ 



TiS) = 
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r 
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s 
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1 



1 + S 
1-S 



\ 

s 







T^J 



u{5, P) = -PS + ^BiiS^ 



u{S,P) = -PS+^BaaS^ 



uiS) = (fill - Bi2)(5" = 2Bee<5^ 



u{S, P) = -3P{S + S^) + i(2Bii + B33 + 2Bi2 + 4Bi3)S^ 



U{S) = 2544^^ 



u{S) = i(2Sii - 2Si2 + B44 - 4S24)5^ 



TABLE III: Calculated stress-strain coefficients and various polycrystalline shear modulus predictions for the stable phase. All 
quantities are in units of Mbar except that volume has been scaled by the ambient volume Vo=13.518 A'^ 



Stable Volume Pressure 



Bulk 



Single crystal stress-strain coefficients 



Polycrystalline 



phase 


iV/Vo) 


(■P) 


modulus {K) 




S33 


B12 


Bl3 


S44 


S24 


VTS 


Voigt 


Reuss 


bcc 


1 


0.00 


1.82 


2.30 


2.18 


1.52 


1.64 


0.51 


0.18 


0.39 


0.41 


0.35 


bcc 


0.779 


0.73 


4.21 


4.79 


4.44 


3.74 


4.09 


0.87 


0.49 


0.43 


0.59 


0.27 


P 


0.754 


0.87 


4.63 


5.52 


5.00 


3.96 


4.43 


0.75 


0.57 


0.51 


0.67 


0.33 


P 


0.729 


1.03 


5.09 


6.00 


5.63 


4.42 


4.84 


0.77 


0.65 


0.50 


0.70 


0.25 


7 


0.705 


1.20 


5.59 


6.17 


6.53 


5.16 


5.27 


1.29 


0.46 


0.69 


0.83 


0.50 


7 


0.659 


1.62 


6.77 


7.87 


7.89 


6.19 


6.23 


1.61 


0.52 


1.03 


1.14 


0.87 


7 


0.636 


1.87 


7.45 


8.79 


8.58 


6.73 


6.85 


1.81 


0.57 


1.19 


1.31 


1.04 


7 


0.593 


2.45 


9.02 


10.81 


9.85 


8.11 


8.38 


2.21 


0.70 


1.43 


1.60 


1.24 


7 


0.572 


2.78 


10.09 


11.56 


9.72 


9.03 


9.97 


2.29 


0.91 


0.99 


1.43 


0.09 
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0.568 


2.84 


10.15 


11.48 


10.64 


9.07 


9.91 


2.04 


1.19 


0.92 


1.37 


0.56 


bcc 


0.551 


3.15 


10.99 


12.71 


11.84 


9.70 


10.56 


2.37 


1.22 


1.29 


1.68 


0.94 
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